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Two-body spectra of pseudoscalar mesons with an 0(a^)-improved lattice 
action using Wilson fermions* 
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We extend our calculations with the second-order tree-level and tadpole improved next-nearest-neighbor action 
to meson-meson systems. Correlation matrices built from interpolating fields representing two pseudoscalar 
mesons (tt-tt) with relative momenta p are diagonalized, and the mass spectrum is extracted. Link variable 
fuzzing and operator smearing at both sinks and sources is employed. Calculations are presented for two values 
of the hopping parameter. The spectrum is used to discuss the residual interaction in the meson-meson system. 



The emergence of hadronic forces from first 
principles is a fundamental question that poses it- 
self to QCD. Remarkably, only few attempts have 
been made to extract effective interactions, or po- 
tentials, between two composite hadrons from the 
lattice . This task is numerically very challeng- 
ing since the residual interaction between color 
singlet composites is about 10"^ to 10"'^ smaller 
than a typical hadron mass. In this paper we 
describe an attempt to study the residual inter- 
action in a simple meson-meson system from four- 
dimensional lattice QCD. 

The improved gluonic action used in our study 
includes planar six-link plaquettes Urt in addition 
to the elementary plaquettes Upi 
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with the hopping parameter k. Both actions are 
corrected for discretization errors to O(a^) at the 
classical level and contain the tadpole factor uq 

ii- 

Starting from a 7r+-like field where the quarks 
carry the flavors u and d 



The improved Wilson fermionic action involves 
nearest-neighbor and next-nearest-neighbor cou- 
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we construct a pion-pion operator jj] with relative 
lattice momentum p and total momentum P = 



(4) 
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To enhance the overlap of the interpolating field 
with the ground state we employ link variable 
fuzzing and operator smearing In particu- 
lar, fuzzy links are used to construct the smeared 
fermionic operators. Smearing is done at both 
sinks and sources. Quark propagator matrix ele- 
ments are computed using a random-source tech- 
nique Q. 

The meson four-point time correlation matrix 

Cf,<t, to) = <i>,-(to)) - {^lit)){^.ito)) (5) 

consists of a free and an interaction contribution 

C^C + Ci. (6) 

The free two-pion correlator is diagonal 

Cpqit, h) = Spq' \cp{t, io)P , (7) 

where is the single-meson two-point function. 
In the present study we are interested in the 
partial wave £ = contained in the irreducible 
representation Ai of the lattice symmetry group 
0{3,Z). We will identify the corresponding re- 
duced matrix elements by a superscript (Ai). 

The eigenvalues of the correlators decrease ex- 
ponentially with increasing time t 
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The mass spectra were obtained from linear fits to 
the logarithm of the eigenvalues of the correlators. 
We define an effective interaction through B 



ot 
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This yields the matrix elements {'p\Hi\q). The 
Fourier transform to coordinate space 
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ip-q\Hi\p + q) 
= ^100(7^) + nonlocal part (10) 
contains a local potential 
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which stems from the p independent part of the 
Hi matrix element in (|l0|). The s-wave projection 
{£ = 0) is 
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where jo is a spherical Bessel function. 

This preliminary study of the two-pion system 



was performed on an L x T — 9 



13 lat- 



tice with /3 = 6.2 in the quenched approxima- 
tion. From the string tension the lattice con- 
stant is determined to be a « 0.4fm correspond- 
ing to w 500MeV. Measurements were taken 
on 48 independent gauge field configurations sep- 
arated by at least 1024 sweeps. Figures 1 (a) 
and (b) show the energy spectra for the free and 
interacting two-pion system according to (||) for 



K — 6.44 and k 
the critical value is k, 



6.57, respectively. Since 
5.5 these translate 



c 

into rather large quark masses. The spectra re- 
veal that the residual interaction lowers the large- 
momentum levels of the free system, thus indicat- 
ing attraction at short relative distances. 

Figures 1 (c) and (d) display the s-wave pro- 
jected local potential resulting from (12). We 
obtain an attractive pion-pion potential with a 
range of < 2a and a depth of about la~^. A 
related study with staggered fermions seems 
to give less attraction. The highest momentum 
p = ^n, n = 4 caused some numerical instabil- 
ity (though no qualitative change) to V\oc{r). We 
therefore only used n = ... 3 in our calculation 
of the effective local potential. The solid lines in 
Figs. 1 (c) and (d) are the results using 48 gauge 
field configurations. The dotted lines represent 
the errors from a jackknife analysis with 8 sam- 
ples of 42 configurations each. The large errors 
in (c) indicate that the results for = 6.44 
are still numerically unstable. For the current 
exploratory analysis no attempts were made to 
estimate systematic errors. 

We expect that a refined analysis, which would 
include extrapolation to the chiral limit and 
proper isospin of the tt-tt interpolating fields, will 
make comparison to experimental findings possi- 
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Figure 1. Energy levels for the free and interacting pion-pion system for k ^ = 6.44 (a) and k ^ = 6.57 
(b) and the corresponding s-wave projected effective local interaction potentials (c) and (d). 



ble. Studies similar in spirit to the one initiated 
here may be extended to various other hadron- 
hadron systems, like tt-N for example. 
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